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Abstract. A S-dimensional quasiperiodic lattice, with overlapping unit cells and 
periodic in one direction, is constructed using grid and projection methods pioneered 
by de Bruijn. Each unit cell consists of 26 points, of which 22 are the vertices of a 
convex polytope 7^, and 4 are interior points also shared with other neighboring unit 
cells. Using Kronecker's theorem the frequencies of all possible types of overlapping 
are found. 



PACS numbers: 05.50.-fq, 61.44.Br 



1. Introduction 

Since the startling discovery of five-fold symmetry in quasiperiodic materials in 1984 
PP, a great deal of research has been done on this subject by both physicists 
and mathematicians. Originally, quasicrystals were constructed by filling the space 
aperiodically with nonoverlapping tiles, such as in Penrose tilings [21IS1I1]|- However, 
recently, Gummelt jB] motivated by physical considerations, proposed a description of 
quasicrystals in terms of overlappings of decagons. Further research [3 El El UDl UH 1121 
IT^ has shown that this may be a more sensible way to understand quasicrystalline 
materials — made of overlapping unit cells sharing atoms of nearby neighbors jH|. Here 
we use a multigrid method to produce a new example of 3-dimensional overlapping unit 
cells. 



2. Multigrid construction 

It is well-known that a Penrose tiling can be obtained by the projection of a slab of the 
5-d euclidean lattice onto a particular 2-d plane V [H [El HE] , and hence its diffraction 
pattern |T71 ^1 has ten-fold symmetry. It is also known that not all lattice points 
k in 7} are allowed (in the sense that they can be mapped onto vertices of a Penrose 

I For a recent review of the theory of quasicrystals, see, e.g., [1]. 
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tiling); only those points whose projections into the 3- dimensional orthogonal space W 
are inside the window of acceptance ^3 120] contribute. The window has been shown 
[TH] to be the projection of the 5-d unit cell Cu(5) with 2^ vertices into this 3-d space 
W. Each facet shared by two neighboring 5-d unit cell cubes is 4- dimensional and 
when projected into 3-d space it produces a polyhedron /C with 12 faces. Therefore, the 
projections of two adjacent 5-d unit cells into 3d must share a common projected facet 
/C. Thus the idea of overlapping decagons must have its extension to three dimension, 
by projecting the 5-d lattice into the space W. 

If dj are the generators of the plane V and Wj are the generators of its orthogonal 
space W, then the projection operators are the matrices = (cZq, . . . , (1^) and 

= {wq, . . . , W4), such that D^W = W^D = 0, where the superscript T denotes 
matrix transpose. More specifically, we choose 

dJ = {cos j9, smj0), wJ = (cos2j0, sm2je, 1) = (dJ., 1), (1) 

where j = 0, ... ,4 and 9 = 2n/5. Using notations and ideas introduced by de Bruijn 
we consider the five grids consisting of bundles of equidistant planes defined by 

X cos 2j9 + y sin 2j9 + z + 'jj = wJR + 7^ = kj, (2) 

for j = 0, . . . , 4, G Z. In (12)), = {x,y,z), and the 7^ are real numbers which 
shift the grids from the origin. We denote their sum by 

7o + 7i + 72 + 73 + 74 = c. (3) 

Without loss of generality, we may restrict c to < c < 1. 

It has been shown by de Bruijn jlj that the Penrose tiling associated with a 2-d 
pentagrid has simple matching rules only for c = 0. For < c < 1 the corresponding 
generalized Penrose tilings do not satisfy simple matching rules, and have seven different 
sets of vertices corresponding to the different intervals of c fIT\ I22j. Nevertheless, the 
diffraction patterns are believed to be the same for all values of c 123 12S] • 

Let the integer kj be assigned to all points sandwiched between the grid planes 
defined by kj — 1 and kj. Then, five integers 

K,{R)=\wjR + j,], j = 0,...,4, (4) 

with [x] the smallest integer greater than or equal to x, are uniquely assigned to every 
point in M'^. A mesh in is now an interior volume, enclosed by grid planes, 
containing points with the same five integers. One next maps each mesh to a vertex in 
Wby 

4 

g{R) = J2k,{R)w, = W^K{R), K^{R) = {Ko{R) , . . . , K^iR)) . (5) 

j=0 

The resulting collection of vertices C = {g{R)\R E M^} is a 3-dimensional aperiodic 
lattice. It has been proven by de Bruijn [T3|, that a point k in satisfies the so-called 
mesh condition and therefore can be mapped into C if and only if D^{k — 7) = D^X, 
where 7""" = (70, . . . , 74) and A""" = (Aq, . . . , A4) with < Aj < 1 so that A is a point 
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inside the 5-d unit cube Cu(5). Thus, the window of acceptance is the interior of the 
convex hull of the points D^rii, where Uj are the 2^ vertices of the 5-d unit cube Cu(5), 
see figure 1. We choose the 32 n^'s as follows 

no = 0, nsi = Co + ei + 62 + 63 + 64, 

IT'j+l = ^l-2j, = ^l-2j + ^4:-2j, IT'j+ll = ^l-2j + ^2~2j , 

= 6l-2j + ^2-2j + ^5-2j, '^'3+21 = ^l-2j + ^2-2j + ^4:-2ji 

nj+26 = ei_2j + e2-2j + 64-2^ + 65-2^, (j = 0, . . . , 4), (6) 
where eg, . . . , 64 are the standard unit vectors in R^, with subscripts counted mod 5 



13 18 




12 19 



Figure 1. The projection of the 5-d unit cube Cu(5) into the orthogonal 2-d space 
T). The window of acceptance is the interior of outer decagon Q whose vertices are 
given by lO- The innermost decagon is denoted by Q and the middle decagon by Q 
with vertices given by There are ten triangles of type (a), which are all further 

subdivided in eight regions of type (al), . . . , (a8), as is indicated for one case in the 
magnification on the left. This is determined by possible overlaps of this triangle with 
the ten regions each of types (b) and (c), see text. 

The projection of Cu(5) with these 32 points into W is a polytope V with 40 edges 
connecting the 22 vertices, and with 20 faces. We let Pi = W^rii for i = 0, . . . ,31. 
The bottom is Pq = (0,0,0) and top is P31 = (0,0,5); they are called the tips of the 
polytope. The remaining twenty vertices of V are 

Pj+l = (dj, 1); Pj+6 = {dj + dj+i, 2), 

P,+2i = {-d,-2 - 3), P,+26 = {-d,^i. 4), (7) 
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for j = 0, . . . , 4. The other 10 points Pn, . . . , P20 are in the interior of the polytope 
and are given by 

Pu+j = {dj + dj^2, 2), Pi6+j = {—dj+i — dj^i, 3), (8) 

again for j = 0, . . . , 4. 

The orthogonal projection of Cu(5) into P is a decagon Q with 10 edges connecting 
the 10 vertices, see figure 1. Let = D^rii for i = 0, ... ,31. Then the vertices of the 
decagon are 

Qn+j = -Pd3-2j, Qie+j = pd5^2j, (9) 
with j = 0, ... ,4, and p = + l)/2. The remaining 22 points Qq, . . . , Q^g 
Q21) • • • ) Q31 in the interior; they are given by 

Qo = Q3I = 0) Qj+l = ^5-2j, Q26+j = -d2-2j, 

Qj+6=P~^d4-2j, Q21+j = -P~^d3^2j- (10) 

Thus if the orthogonal projection D^(k — 7) is in Q, then its projection W'^k is in C 
Consider a polytope P in £ whose bottom is the projection of a point k which is 
orthogonally projected onto the center of the decagon, then the ten vertices of V at 
height z = 1 or z = 4 correspond to the vertices of the middle decagon Q, see figure 1; 
the ten vertices of the polytope at height z = 2 or z = 3, correspond to the vertices of 
the innermost decagon Q; while the ten interior points of the polytope correspond to 
the vertices of the (outer) decagon Q shown in figure 1. 

The points in the quasiperiodic lattice C need careful analysis. We shall show that 
C consists of polytopes V with 22 vertices and 4 interior points which are shared with 
neighbouring polytopes. We shall also show that it is periodic in the z-direction — the 
direction of the line joining the two tips of V — with period 5 corresponding to the height 
of V. 

3. Overlapping unit cells 

The window Q is known |31 EE] to be everywhere dense and uniformly distributed. Thus 
each point in Q corresponds to a point in the quasiperiodic lattice C Using an idea of 
de Bruijn jl], we may find out the condition for both k and fc + to satisfy the mesh 
condition — to lie both in the window of acceptance, which give insights to how points 
in C are related to one another. 

It is straightforward to show that every point inside the innermost decagon Q 
corresponds to a point in C that is connected with 10 neighbors, and is in fact a tip of a 
polytope, as the middle decagon Q whose center is shifted to a point inside Q still lies 
inside Q. This innermost decagon Q is further divided into 10 triangles §. Whenever 
the center of a decagon Q is shifted to a point inside one of the triangles, four of the 

§ A rotation of points in Q by an angle of 2tt£/5, corresponds to a rotation of the polytopes in C by 
an angle of -iTTi/b about the z-axis, while a rotation by tt — inversion of points in Q through the origin 
of decagon Q — corresponds to the 3-d inversion of polytopes in C through the center of each polytope. 
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vertices of the shifted outer decagon now he inside Q. More precisely, if D^{k — 7) is 
inside triangle (a) in figure 1, then W^k is a tip of a polytope V whose interior points 
W^{k + n2o), W^{k + rirs), W^{k + riis) and W^{k + nn) — corresponding to the 
four points on Q on the opposite side of triangle (a) — are now also in C. This means 
each polytope in £ can have only four interior points which are also in £. Thus each 
such unit cell contains 26 atoms, 22 exterior, and 4 interior sites. 

It is also easy to find out how the polytopes share these interior points. This 
is equivalent to finding the condition that both k and k + Ui ior i — 1, . . . , 10 or 
i = 21, . . . , 30 are orthogonally projected into Q. 

Consider the ten vertices of the polytope at height z = 1 or z = 4 corresponding 
to the vertices of middle decagon Q. When the center of the decagon is shifted to a 
point in one of ten rhombs of type (b) shown in figure 1, then one of the ten vertices 
of the shifted Q is inside Q. For the example in figure 1, VF^(fc) and VF^(fc + n28) 
are both tips of polytopes, with the point on Q on the opposite side of rhomb (b). 
As a consequence, the polytope whose top is W^(fe + n28) shares with the polytope V 
whose bottom is W^(fe) a polyhedron /C with 12 faces. 

Moreover, the vertices of the polytope at 2; = 2 or 2; = 3 correspond to the vertices 
of Q. When the center of this decagon is shifted to a point in one of ten octagons of type 
(c) shown also in figure 1, then one of the ten vertices of the shifted Q moves inside Q. 
For the example in figure 1, W'^{k) and W'^{k + 7225) are then both tips of polytopes, 
while 71,25 is the point on Q on the side opposite to the octagon. Such two polytopes 
share a polyhedron J' with 6 faces. When vertices of V at both z = 3 and z = 4 are 
inside Q, they are tops of some polytopes|| that share with V either a polyhedron JC or 
J. 

By considering how these twenty regions intersect, we find that each of the triangles 
in Q is further divided into eight regions, as shown for one case in figure 1, with this 
triangle magnified on the left. When the orthogonal projection of fc — 7 is in (al), 
its projection into W is a polytope V intersecting with four others whose tips are at 
W'^(k + 7x21), W^{k + Us), W^{k + 7224) and W'^{k + tiq), and sharing with each of 
them a polyhedron of type J'. When the projection is in (a2) or (a3), either {k-\-n26) 
or + 7ii) becomes also a tip of a polytope, so that V intersects with five polytopes, 

sharing with one of them a polyhedron of type /C. In (a8), V intersects with all of the 
above six polytopes. In regions (a5) or (a7), (fc + 7i2i) or W'^{k + nQ) is no longer a 
tip and V intersects with five polytopes sharing with two of them a polyhedron of type 
JC. Finally, while in (a4) or (a6), either W'^{k + rii) or W'^{k + n26) is no longer a 
tip, and V intersects with four polytopes sharing with one of them a polyhedron /C and 
with the other three a J^. 

Therefore, in summary, there are only five different possibilities: 

(i) The points inside a quadrilateral of type (al) correspond to a polytope intersecting 
with four other polytopes sharing with each a polyhedron of type J^. An example 

II We shall later show that lattice C is periodic in the ,2-direction with period 5. 
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of this case is shown in figure 2a. 

(ii) Points inside triangles of type (a2) or (a3) correspond to a polytope intersecting 
with five other polytopes, sharing with one of them a polyhedron K, and with the 
other four polyhedra of type J . Such a case is shown in figure 2b. 

(iii) If the point is inside triangles of type (a4) or (a6), the polytope intersects with four 
other polytopes sharing with one of them a polyhedron /C and with the other three 
polyhedra of type J . 

(iv) If the point is in a triangle (a5) or (a7), the polytope intersects with five other 
polytopes sharing with two of them polyhedra of type K, and with the other three 
polyhedra J . 

(v) Finally, if the point is inside a pentagon (a8), the polytope intersects with six other 
polytopes sharing with two of them a /C and with the other four a J , as is shown 
in figure 3. 

The relative frequencies are related to the ratios of their areas and therefore the 
normalized probabilities are given by 

Pal = 2p-^ Pa2 = Pa3 = P~\ P.A = P.6 = P"', 

Pa5 = P.7 = P'\ P.S = + P-'. (11) 



4. Parallelepipeds and windows 

Now the condition that k satisfies the mesh condition can be further simplified by 
considering parallelepiped P^k^, k^, ki) sandwiched between the six grid planes k^ — 1, 
k^, ko — 1, ko, ki — 1, and ki. Using Q for the grid planes, we can solve for the points 
of intersections, and find that for every point R in P{k4^, k^, ki), we may write 

Ko{R) = ko, Ki{R) = ki, K2{R) = [a\ +ki + m, 

K^{R) = [(3\ + A;i + n, K^{R) = k^, (12) 

where \_x\ is the greatest integer less or equal to x, while m and n are integers satisfying 
— 1 < m, n < 2, and 

a = p^^{ko - /ci - 7o + 7i) + 72 - 74, 

/3 = p-i(A;o - - 7o + 74) + 73 - 7;^. (13) 

The range of m and n in equation (fT^ is limited by their possible values at the eight 
corners of the parallelepiped, but not all 16 choices are allowed by the mesh condition. 
Therefore, we now project K{R) given by (|T2|) to V and find 

4 

D^{K{R) - 7) = J2{K,{R) - 7,)d, = (m - a)d2 + {n - b)d,, (14) 

j=0 

in which a = {a} = a — [aj and b = {(3} = f3— lf3\ . The vector ()14|) must lie within the 
decagon Q. Thus the allowed values of m and n are determined by a and b only. As the 




(a) (b) 

Figure 2. (a) When the orthogonal projection of fc — 7 is in region (al), its projection 
in £ is a polytope sharing an interior point with each of four neighboring polytopes 
whose tips are projections of fe + n2i, k + ns, k + 9124 and k + n^; the shared interior 
points are projections of k + rin, fc + nig, k + ni4 and fe + nig. (b) When it is in region 
(a2), the polytope has five neighboring polytopes sharing the same interior points. The 
additional polytope has its top at {k + n26)- 



differences ko — and k^ — ki run through all integer values, we find from Kronecker's 
theorem j2Z| that a and b are everywhere dense and uniformly distributed in the interval 
(0, 1). Proofs of such "ergodicity" in more general situations can be found, e.g., in the 
works of Hof and Schlottmann |2H] • 

5. Periodicity in third direction 

Furthermore, because of the difference property shown in ()13|). a and (3, which 
determine the configuration of the parallelepiped, remain the same, if (^0,^1,^4) —>■ 
{ko + i,ki + i,k4 + i). As a consequence we find K{R) — > K{R) + in^i, and its 
projection is periodic in the z direction with period equal to 5. It is also interesting to 
note that if 7 ^ 7 — |cn3i, a and P in (jlHj) are also unchanged, so that the projection 
into the 3d space does not show drastic changes when c 7^ 0, which is behavior very 
different from the 2-d case 21j. 

It is not difficult to find values of (a — m, 6 — n) in (fT^ corresponding to the ten 
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(a) (b) 

Figure 3. The polytope V and its six neighbors in C for case (a8). Their intersections 
with V are indicated with different colors. In (a) the orientation is (0, 0) — (0, 80°); in 
(b) {6,(1)) = (90°, 0), so that each polytope appears as a decagon. 

vertices of the decagon Q given by 0- We find 



Ql6 ^ 


(a-0,6- 


0) = 


(1,1), 


Qii ^ 


(a-0,6+ 1) 


= {0,p), 


Ql8 


(a - 1,6 + 


1) = 


(-1,P), 


Ql3 ^ 


(a -2, 6+1) 


= 1) 


Q20 ^ 


(a-2,6- 


0) = 


(-p,o). 


Ql5 ^ 


(a - 1,6- 1) 


= (0,-1) 


Ql7 ^ 


(a — 1, 6 — 


2) = 


(0,-p), 


Q12 ^ 


(a -0,6 -2) 


= (1,-P) 


Q19 ^ 


(a + 1,6 - 


2) = 


(P,-l), 


Ql4 ^ 


(a +1,6-1) 


= (p,o). 



up to ambiguities when a or 6 is integer, as m or n changes by 1 when choosing a or 6 
to be or 1. The edges of decagon Q lead to hnear equations in {a — m,b — n) and the 
mesh condition for K{R) becomes a set of inequahties in a — m and 6 — n as shown in 
figure 4. Hence, it is very easy to create a routine to generate C 

We note that C is periodic in the longitudinal or z-direction with period 5, and 
aperiodic in the horizontal directions. This is a model with five layers which repeat in 
the z-direction periodically. In each of the layers, the points behave similarly, and the 
allowed sites for 2; = 3 are shown in figure 5 (a). This shows decagonal symmetry, but 
the layer may not be covered by Penrose tiles. In figure 5 (b), we have sites for atoms 
in all five layers plotted, with sites in different layers represented by different symbols 
and colors. 

Therefore, C is another example of a decagonal quasicrystal. Decagonal 
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Figure 4. The image of Q under the inverse mapping of H14|l into the {a — m,b — n)- 
plane. The boundaries of regions with different (m, n) are given by linear equations in 
terms oi x = a — m and y = b — n. 
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(a) (b) 

Figure 5. The allowed positions in each layer, (a) Allowed sites at height z = 3 are 
plotted, (b) All layers are plotted: points for z = 1 are represented by black boxes; for 
z = 2 by red circles; for z — 3 by blue diamonds; for z = 4 by green crosses; for z — 5 
by brown points. 



quasicrystals exist in nature, and have been extensively studied experimentally and 
theoretically [201 EH E21 ESI EH EH] • From Steurer's review [HI], one finds translational 
periods along the tenfold z-axis varying from 4 A to 16 A for different alloys, which may 
allow 2 layers [^21; or even 5 or more layers, within a period. 

A situation somewhat similar to ours is found by Ben-Abraham, Lerer and Snapir 
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[T^ . who find that the projection of a 6-d lattice to 3d for a certain choice of bases 
produces a quasicrystal lattice which is also periodic in the z-direction. Their model 
has six-fold rather than five-fold symmetry. 

The projection of — is well-known |3m |H3 EHl- However, interpreting the 
3-d quasicrystal lattice C as overlapping polytopes V provides a more systematic way to 
understand this lattice. We find that each unit cell has 26 sites, sharing the four interior 
sites with its neighbors. The lattice is periodic in the z-direction with period 5, and 
quasiperiodic in the xy-directions. Unlike the projection of — > M^, which undergoes 
drastic change in behavior (no inflation and deflation rules) if c 7^ 0, the projections 
into 3d are in the same class for c = and c 7^ 0. 

Since real quasicrystals have icosahedra or triacontahedra as unit cells, which are 
projections of a six-dimensional hypercube to a 3-d space [101, the above method perhaps 
can also be used to determine all possible overlappings and their frequencies. 
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